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Abstract — The paper consider a complex dynam- 
ics of electron beam with virtual cathode and local 
neutralization of the beam charge density near an- 
ode. Different types of nonlinear behaviour, includ- 
ing deterministic chaos, were treated. It is shown 
that chaotic dynamics arises as results of spatiotem- 
poral structures interaction. 

Introduction and model 

In the present paper we consider a electrostatic short- 
circuited diode model. The immobile ion background 
with concentration n p locates near left boundary (an- 
ode plasma). The electron beam with overcritical cur- 
rent injected into diode with nonpurtubated velocity 
vq and charge density po. ln this case virtual cathode 
(VC) forms in the beam as a result of the electrostatic 
instability 0, and some beam part is reflected from 
VC to injected boundary. VC oscillates in the diode 
region. The electron beam with VC demonstrates wide 
diversity of nonlinear phenomenons, including chaotic 
behaviour 0, synchronization || and other. Investi- 
gation of complex dynamics attracts many researchers, 
that is, such behaviour is the characteristic property of 
the beam with VC. That counts is study the structure 
formation, since it is well known, that chaotic dynam- 
ics in distributed systems is connected with pattern 
formation (see 0] and references therein). 

The such simplest model of device with VC as a pla- 
nar diode with overcritical current describes a different 
nonlinear phenomenons in the real vircator systems. 
Our model with local neutralisation is a simple model 
of vircator with injected plasma 

The behaviour of system is determined by the di- 
mensionless parameter related current 

a = uj p L/vq, 

where oj p is the beam plasma frequency, L is the dis- 
tance between diode planes and neutralisation param- 
eter 

n = Up/no. 



Hence no = po/e and value of plasma region length x p 
is constant (x p = 0.25L). 

The effect of neutralisation degree of anode plasma 
on VC dynamics was investigated with the aid of 
particle-in-cell simulation. The macroparticles in the 
simulation obey the non-relativistic equations of mo- 
tion 

dx/dt = v, 

dv/dt = —(q/m)d(f)/dx 1 

where x is the position of particles, v is the velocity of 
of the particles, q is the charge and m is the mass of 
the macroparticles. The code integrates the equation of 
motion forward in time using a leapfrog scheme. The 
potential 4> is computed by the Poisson's equation in 
one dimension 



d 2 <j>(x) 
dx 



-a 2 {p{x) - p p (x)). 



Hence p{x) is the spatial distribution of beam charge 
density and p p {x) is the distribution of immobile ion 
background. In our case 



p p (x) 



X < X 



pi 



System dynamics 

The tentative analysis of nonlinear dynamics were effec- 
tuated from observation of time series of electric field 
oscillation in the injection plane. Power spectra and 
projections of attractors were reconstructed from time 
series. Based on this analysis, domains for distinct be- 
haviour were isolated in parameter plane (a, n) (see 
Fig. 1). 

The VC oscillation (VCO) for small value of neutral- 
isation is regular (domain marked A in the parameter 
plane; Fig. 2(a)). Analysis of physical processes shows 
that only one electron bunch (VC) is arisen in the sys- 
tem. This bunch is marked on the spatiotemporal di- 
agram (Fig. 3(a)). Besides, metastable particles, which 
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Fig. 1: Bifurcation diagram on the parameter space 
(a, n). White area S corresponds to nonuniform equi- 
librium 



exist in the interaction space during of more than one 
period of VCO, is observed in the beam. However, 
charge density of the metastable bunch is small, and 
it is little affected by VCO. The weakly chaotic VCO 
arises as neutralisation and current increase (domain B 
for n < 2.0; Fig. 2(b)). In this case metastable bunch 
density grows. The further increasing of n leads to 
formation of profound metastable bunch in the beam 
(Fig. 3(b)). A buildup of space charge density in VC 
region entails the regime with large base frequency in 
the VCO spectrum (compare Fig. 2(a,b) and Fig. 2(d), 
that obtain for the same value of a = 2.12571"). This 
behaviour of system take place for n > 2.25 -r 2.5 (B; 
Fig. 2(d)). A change-over from the weak chaos for 
small values of n to the weak chaos for large n de- 
rives through two domains of strongly chaotic VCO. In 
the first regime (domain C in Fig. 1) phase portrait is 
homogeneous, there are not sharp peaks in the power 
spectrum. The second regime (D) may be treated as 
intermittence (Fig. 2(c)). 

For large values of n (domain E) system demon- 
strates highly non-regular oscillation with noise-like 
spectrum and homogeneous attractors (Fig. 2(e)). In 
this case VC is formed out of anode plasma region. VC 
exists constantly and chaotic dynamics is determined 
by the reflection of the particles from VC. Note, VC is 
not moved in the space, but depth of potential barrier 
is oscillated in time. 

Dimensions of the reconstructed attractors was esti- 
mated for different types of chaotic behaviour. Fig. 4 
presents correlation dimension of attractors D versus 
value of embedded dimension m for strong chaos. Di- 
mension is saturated for small values of m. Small values 
of embedded dimensions justify appearance of chaotic 
behaviour in the beam in the result of interaction be- 
tween small numbers of structures. 
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Fig. 2: Power spectra, reconstructed phase portraits 
and time series for different regimes 
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Fig. 3: Spacetime diagrams for regular (a) and chaotic 
(b) oscillation 
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Fig. 4: Dependence of attractor dimension from em- 
bedded dimension 



Pattern formation 

The spatio-temporal data of charge density p(x, t) were 
analyzed by the Karhuncn-Locvc orthogonal decom- 
position |||. This method decomposes a data set into 
spatial orthogonal modes {ip(x)}i. This modes is the 
solution of the integral equation 



R(x,x')i>(x')dx' = Ai/j(x) 



where R(x,x') = {p(x , t) p(x' , t)) t is the mutual corre- 
lation function. The value of eigenvalue is propor- 
tional to the energy of i-th mode. Karhuncn-Loeve 
method is optimal in the sense that one is optimized 
(from the viewpoint of energy of modes) eigensets {A}.; 
and The measure of energy of modes is 

Wt " 5V 

The energy of several first modes for different values 
of n is presented in Table for a — 1.357T. The typi- 
cal spatial distributions of modes are shown in Fig. 5. 
In the all regimes of VCO behaviour of beams is de- 
termined by the small number of structures, because 
more then 90% energy is contained in the 3-^4 higher 
modes. 

Table: Energy Wi (in %) of Karhunen-Loeve modes 



Number of 


Value of neutralisation, 


n 




mode, i 


0.0 


0.35 


0.7 


1.05 


1.4 




1 


67.1 


62.1 


60.0 


56.1 


52.6 




2 


17.1 


16.4 


20.0 


21.8 


23.3 




3 


5.7 


8.3 


7.1 


7.9 


8.5 




4 


3.1 


4.6 


4.1 


4.6 


5.0 






















Value of neutralisation, n 






1.75 


2.1 


2.45 


2.8 


3.15 


3.5 




52.3 


60.8 


71.5 


70.3 


69.8 


83.2 




25.1 


21.4 


10.2 


12.3 


12.7 


11.5 




8.5 


7.4 


6.3 


4.2 


4.1 


3.0 




4.6 


2.9 


2.2 


3.3 


3.6 


0.4 



For small neutralisation first modes demonstrates 
strong nonuniform charge density distribution with one 
peak (Fig. 5(a)), that corresponds of typical distribu- 
tion of density in VC. Second mode describes are pro- 
cesses of destruction of VC and disposal of particles 
from VC to anode. Third and other modes correspond 
to additional bunches in the beam. Increasing of n 
leads to growth of second modes at the expense of 
higher mode. 

Metastable bunch is formed for values of n ~ 2 3. 
In this case first and second modes together describes 
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Fig. 5: Typical structures for different values of neu- 
tralisation 



dynamics of VC and metastable bunch, besides cross- 
correlation between dynamics modes A\(t) and ^b(i) 
is large. Hence temporal dynamics of modes is 

Ai(t) = J p(x,i)tpi(x)dx. 

Spatial distribution of modes is strongly localizing 
in the interaction space for regime E (Fig. 5(c)). Only 
one structure — weakly oscillating in the space virtual 
cathode — exists in the beam, and more than 80% 
energy accumulated in the first mode. Chaotic VCO is 
determined by the reflection of the larger part of beam 
from VC. 

Conclusions 

In the electron beam with VC and local neutralisation 
different types of nonlinear oscillation are recognized. 
Influence of density of anode plasma on the chaotic 
dynamics of VC are considered. Large neutralisation 
degree leads to strong chaos in the VCO. Strange at- 
tractor is most homogeneous in this case. Relationship 
between different type of chaotic behaviour and struc- 
ture formation are shown with the help of analysis of 
physical processes in the diode. The typical patterns 
(VC and different additional structures) were recog- 
nized for regimes with small and large neutralisation. 
It is shown that appearance of chaotic behaviour con- 
nects with growth of charge density in the additional 



structures. Investigation of chaotic behaviour in our 
model may be helpful for design of vircators with an- 
ode plasma grid, because changing of density of anode 
plasma lead to changing of radiation characteristics in 
these devices. 
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